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Abstract. In this paper, we introduce a class of new generalized super Bell 
polynomials on a superspace, explore their properties, and show that they are 
a natural and effective tool to systematically investigate integrability of super- 
symmetric equations. The connections between the super Bell polynomials and 
super bilinear representation, bilinear Backlund transformation, Lax pair and 
infinite conservation laws of supersymmetric equations are established. We take 
supersymmetric KdV equation and supersymmetric sine-Gordon equation to il- 
lustrate this procedure. 
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1. Introduction 



The supersymmetry represents a kind of symmetrical characteristic between 

boson and fermion in physics. The concept of supersymmetry was originally 

introduced and developed for applications in elementary particle physics thirty 

years ago [T] [3] . It is found that supersymmetry can be applied to a variety 
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of problems such as relativistic, non-relativistic physics and nuclear physics. In 
recent years, supersymmetry has been a subject of considerable interest both in 
physics and mathematics. The mathematical formulation of the supersymmetry 
is based on the introduction of Grassmann variables along with the standard ones 
[36] . In a such way, a number of well known mathematical physical equations have 
been generalized into the supersymmetric analogues, such as supersymmetric 
versions of sine-Gordon, KdV, KP hierarchy, Boussinesq, MKdV etc. It has 
been shown that these supersymmetric integrable systems possess bi-Hamiltonian 
structure, Painleve property, infinite many symmetries, Darboux transformation, 
Backhand transformation, bilinear form, super soliton solutions and super quasi- 
periodic solutions [I] -[IB]- In our present paper, we investigate the integrability 
of supersymmetric equations by using a class of super Bell polynomials which are 
a multidimensional and super generalization of ordinary Bell polynomials. 

The ordinary Bell polynomials introduced by Bell during the early 1930s are 
a class of exponential polynomials, which are specified by a generating function 
and exhibit important properties [22]. The Bell polynomials have been exploited 
in combinatorics, statistics and other fields [23]- [25]. Some generalized forms of 
Bell polynomials already appeared in literature [26]- [30]. More recently Lam- 
bert, Gilson et al found that the Bell polynomials also play important role in 
the characterization of bilinearizable equations. They presented an alternative 
procedure based on the use of the properties of Bell polynomials to obtain pa- 
rameter families of bilinear Backlund transformation for soliton equations. As a 
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consequence bilinear Backhand transformation with single field can be linearize 
into corresponding Lax pairs [3T]-|33j. 

Our paper is a further contribution to the theory of Bell polynomials and 
supersymmetric equations. We reconsider Bell polynomials in a more extended 
context- superspace. We define a kind of new generalized super Bell polynomials 
and discuss their relations with super bilinear equations, which actually provides 
an approach to systematically investigate complete integrability of supersym- 
metric systems. As illustrative examples, the bilinear representations, bilinear 
Backhand transformations, Lax pais and infinite conservation laws of the super- 
symmetric KdV equation and supersymmetric sine-Gordon equation are obtained 
in a quick and natural manner. 

The layout of this paper is as follows. In Section 2, we briefly recall elemen- 
tary notations about super differential, integrals and super bilinear operators on 
superspace. As In Section 3, we propose theory of super Bell polynomials and 
establish their connections with supersymmetric equations. As consequence a ap- 
proach to investigate integrability of supersymmetric equations is presented. In 
the Sections 4 and 5, as applications of super Bell polynomials, we study integra- 
bility of supersymmetric KdV equation supersymmetric sine-Gordon equation, 
respectively. At last, we briefly discuss further possible generalization and appli- 
cations of Bell polynomials and future work in Section 6. 

2. Derivatives and bilinear operators on superspace 
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To make our presentation easily understanding and self-contained, in this 
section we first briefly review some notations about superanalysis [31]- [37] and 
super-Hirota bilinear operators [T9l [20] . 

A superalgebra is a Z 2 -graded space A = A © Ai in which, A is a subspace 
consisting of even elements and Ai is a subspace consisting of odd elements. A 
parity function is introduced for homogeneous elements on the A, namely, \a\ = 
if a G A and \a\ = 1 if a G Ai. 

The superalgebra is said to be commutative if the supercommutator [a, b] = 
ab — (— l)' a " 6 '6a = 0, for arbitrary homogeneous elements a, b G A. 

A commutative superalgebra A with unit e = 1 is called a finite-dimensional 
Grassmann algebra if it contains a system of anticommuting generators 9j,j = 
1, • • • ,n with the anticommutative property: [9j,9k] = OjOk + 9^9 j = 0, 0] = 

0, j,k = 1,2, ••• ,n. 

Let A = A ©A! be a finite-dimensional Grassmann algebra, then the Banach 
space M^' n = A™ x A™ is called a superspace of dimension (to, n) over A. In 
particular, if A = C and Ai = 0, then IR™' n = C m . We may take even- valued 
complex space A = C in our context. 

A function f(x, 6) : M™' n — > A is said to be superdifferentiable at the point 
(aj, 6) G M™' n with even coordinates ) and odd coordinates 6 = 

(#!,■■• ,9 n ), if there exist elements Fj(x, 6), F^x, 6) G A, j — 1, • • • , to; k — 

1, • • ■ , n, such that 

m n 

f{x + h,0 + h) = f(x, 0) + ^i) + + °d K h > ^)ID» 

j=i fc=i 
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where the vectors h = (hi,-- - ,h m ) G A™ and h = (hi,-- - ,h n ) G A™. The 
Fj(x, 6), F k (x, 6) are called the super partial derivative of / with respect to Xj, 9 k 
at the point (x, 6) and are denoted, respectively, by 

9 -^> = Fj(x, 6), d -^> = F k (x, 0), j = 1, • • • , m; k = 1, • • • , „. 

The derivatives l9 ^' 6> - ) with respect to even variables Xj, j = 1,2, ---m are 
uniquely defined. While the derivatives to odd variables 9 k , k — 1, 2, • • -n 

are not uniquely defined, but with an accuracy to within an addition constant 
c9i ■ ■ • 9 n , c e A from the annihilator ± L n , L n — {9i • • • 9 n , 6 e A"}. 

Let / = f(x,0),g = g(x,8) : IR™' n — > A be a superdifferentiable function, 

then super derivative also satisfies Leibnitz formula 

d Xj (fg) = (d Xj f)g + f(d Xj g), j = 1, • • • ,m, 

d 9k (fg) = (dej)g + (-l)^f(d dk g), k = 1, • • • , n. 
Let differential operators V k = dg k + 9 k d Xr ( k — 1, • • • , n; r e {1, • • • , m} ) 

be supersymmetric covariant derivatives, we can show that they satisfy 

V k (fg) = (V k f)g + (-lf\f(V k g), 

(2.1) 

[V v V k } = 0, V\ = d Xr . 
Denote by ^(A^A) the set of polynomials defined on A" with value in A. 

We say that a super integral is a map I : P(A",A) — > A satisfying the following 

condition is an super Berezin integral about Grassmann variables 

(1) A linearity: I(fif + vg) = fil(f) + ul(g), fi, v e A, /, g e V(k n 1 , A); 

(2) translation invariance: I(f^) = 1(f), where f^ = f(0 + £) for all ^ G A", 

feV(k n i,k). 

We denote I(9 £ ) = I £ , where e belongs to the set of multiindices N n = {e = 
(si, •••,£„), Ej = 0, 1, 6 s = 91 1 ■ ■ ■ 9 6 n n ^ 0}. In the case when I e = 0, e G N n , \e\ < 
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n — n — 1, such kind of integral has the form 

/ (/ ) = j( / )/(i,...,i) S ^|-/(i,-,i), 

Since the derivative is defined with an accuracy to with an additive constant form 
the annihilator ± L n , it follows that J : V — > t^j^L n is single- valued mapping. 
This mapping also satisfies the conditions 1 and 2, and therefore we shall call it 
an integral and denote 

J(f) = j f(0)dO = J 9, ■ ■ ■ n d0i • • • d9 n , 



which has properties: 

Q\ • • • Q n dQ\ • • • d6 n = 1, 
df 



dd x ---dd n = 0, j = 1, • • • ,n. (2.2) 

For a pair of Grassmann- valued functions f(x,0), g(x,0) : IR™' n — > A, the 
ordinary Hirota bilinear operator is defined by 

D x .f(x,0) -g{x,0) = (d Xj - d x , j )f(x,O)g(x',0')\ x > =Xt0 > =0 , j = !,■■■ ,m, 

and super-Hirota bilinear operators are defined as 

S k D Xj f(x, 0) ■ g(x, 0) = (V k - V' k )f(x, 0)g(x', 0') \ x , =x , g >=9, k = 1, - ■ ■ ,n, 

here we have denoted V' k = d e > + 9'd x . 

It can be shown that these super-Hirota bilinear operators have properties 
S 2 k N f-g = Dlf-g, NeZ 

S k f-g = (V k f)g-(-lf\f(V k g). 
In our context, we are interested in bosonic, also called even, superfield func- 
tion f(x, 0) : M™' r 1 ->■ WfC = A . It can be expanded in powers of odd coordinates 
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9k, k = 1, ■ ■ • , n, that is, 

k>0 ji< - <j k 

where the coefficients fj 1 -j k (x) G A are even functions with respect to xi, ■ ■ ■ , x m . 

3. Generalized super Bell polynomials on superspace 

Based on the above fundamental notations, in this section we develop theory 
of generalized super Bell polynomials, which are a main tool to study the inte- 
grability of supersymmetric equations. 

2.1. Generalized super Bell polynomials 

To well compare our super Bell polynomials with ordinary ones, let's first sim- 
ply recall the rdinary Bell polynomials. During the early 1930s, Bell introduced 
three kinds of exponential polynomials [22]. 

The first Bell polynomials are defined as 

Ux^r) = e- txr d n x e tx \ (3.1) 

where r > is a constant integer, n > an arbitrary integer, and x,t EM. inde- 
pendent variables. For r = 2, the Bell polynomials £ n (x,t) are exactly Hermite 
polynomials. 

The Bell polynomials are algebraic polynomials in two elements x and t. The 
first few lowest order Bell Polynomials are 

f (x, t, r) = 1, ^(x, t, r) = rtx r ~ l , ^ 2 {x, t, r) = r¥i 2r " 2 + r(r - l)tx r ~ 2 , 

&{x, t, r) = r¥i 3r " 3 + 3r 2 (r - l)t 2 x 2r ~ 3 + r(r - l)(r - 2)tx r ~ 3 . 
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The second Bell polynomials are a generalization of the Bell polynomials (3.1) 
and defined by 

0n = 0(«1, •••,««), 00 = 1, 0n+l = ^ tt s+ i0„_ s , 

where («i, • • • , a n , ■ • • ) is an infinite sequence of independent variables. For the 
particular sequence aj = j\ {^j x t r ~ j \ J = 1, • • • , aj = 0, j > r, we have 

0n = in(x } t,r). 

The Bell polynomial n is a polynomial about variables cm, • • • ,a n . For 
instance, the first three of second Bell polynomials read 

0o = 1, 0i = «i + «2, 03 = a? + 3aia 2 + «3- 

The third Bell polynomials, further generalization of the £ n and n , are de- 
fined by 

Y n = Y n (y t ,--- ,y nt ) =e~ y d?e\ (3.2) 

where y = e at — a = ait + a 2 t 2 /2\ + • • • , and we have denoted derivative notation 
Vkt — dfy. For the spacial case when = 0, k 7^ r, then we have 

Y n £ n (x, t, r). 

The polynomials (3.3) are polynomials about the derivatives of function y, 
for example, the first three are 

Y = 1, Yi = y t , F 2 = 2/t + ?/2t, >3 = Vt + tytVit + Vst, 

More recently Lambert et al generalized the third Bell polynomials as 

Y n = e-*d£---dZ#, (3.3) 

where y = y{x u ■ ■ ■ , x m ) : R m -»■ R |31]-[33]. 
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We now propose the following mult i- dimensional and super extension to the 
ordinary Bell polynomials (3.1)-(3.3). 

Definition 1. Let / = f(x,0) : R™'™ — y A be a superdifferential bosonic 
function, the generalized super Bell polynomials (super F-polynomials) is defined 
as follows 

Y £ . x , e (f) = V<-,,o fr.,,,,.0 = e~fV x ■ ■ ■ V n d x \ • • • , (3.4) 

where £j > 0, j — 1, • • • , m denote arbitrary integers. To make subscript in ex- 
pressions simple, we use some abbreviation notations in our context, for example, 

£. (^1; j @>m)i & (^1) i @n)i £ ' X (£i3Ji, , £m x ra)i 
T ('"lj ' ' ' i ^rn)i T ■ X [t\Xi ) , T m X m ), 

H = (//i, • • • , fl n ), fJ, = (nxBx, • • • , nJn)- 

Remark 1. The first notation Y t . Xt e{f) in (3.4) denotes the ^ -order deriva- 
tives of / with respect to the variable Xj,j = 1, • • • ,m and covariant deriva- 
tives with respect to 9 k , k = 1, • • • , n. The second notation Y^. xQ [f r . Xtfl .g\ 
implies that the super Bell polynomials (3.4) should be understood as such a 
multivariable differential polynomial with respect to partial derivatives f r x ,n-e ( 
Tj = 0, • • • ,£j, j = 1, • • • , m, iik — 0, 1, k — 1, • • • ,n), but not variable elements 
x ji @k (j = I, • • • ,m;k — 1, • • • , n) as ordinary Bell polynomials. For instance, 
the Y 3x (f) in the next example is a polynomial Y 3x (f x , fo x , fox) with respect to 
three variable elements f x , fo x , fox- 

To better understanding our generalized super Bell polynomials, let us see an 
illustrative example. For the special case / = f(x,0i,02), the associated super 
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Bell polynomials denned by (3.4) read 



Y x (f) — fx, Y 2x (f) — hx + fl, 

Y:iAf) f,r ■ -yr,k, • /i 5 . Y dl (f) = v l f, 

YoMf) = ViV 2 f + (VJ)V 2 f, Y xA (f) = V x f x + f x VJ, 
Yzr.olf) = hxVif + VJ 2x + f;V.f + 2f x VJ x , 

Yum if) = fcxVif + Urfi,-P:f + t>f 2x V x f x + 3/^^ + 3f*V x f x + V 1 f 3x . 



Let see the relations between our generalized super polynomials and ordinary 
Bell polynomials, as well as ordinary generalized Bell polynomials. 

For the special case M™' n = M 2 , £ 2 = 0, / = f(xi, x 2 ) = x 2 x\ with the constant 
integer r > 0, then (3.4) reduces to the first Bell polynomials (3.1) 



For the case M.™' n = M. m , the corresponding generalized super Bell polynomials 
(3.4) degenerates to generalized Bell polynomials (3.3) given by Lambert et al. 
The Bell polynomials admit partitional representation [3T] 



where the sum is to taken over all partitions [(r^i, • • • , r ml ) ci , • • • , (r lfc , • • • , r mk ) Ck ] 
the m-tuple (li, - ■ • ,£ m ). 

In following, we investigate properties of super Bell polynomials which are 
key results to establish connects with supersymmetric equations. 

Theorem 1. Under the Hopf-Cole transformation f = In if), the generalized 
super Bell polynomials Y r . XifJ/ .e(f) can be "linearized" into the form 



Ye ixi (f) 



e -x2x{ d ei e 




(3.5) 



(3.6) 
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Proof. According to the definition (3.4), we have 

which finishes the proof of Theorem 1. □ 

Remark 2. According to the theorem, under the Hopf-Cole transformation 
/ = lnip, a equation in term of linear combination of super Bell polynomials, i.e. 

can be linearized into the form 



^2 CrA x > 0)ipr-x,n-e = 0, 



where C r ^(x,6) are functions independent of the function /. This is a key 
property to construct the Lax pair of supersymmetric equations. 

Theorem 2. The super Bell polynomials (3.4) admit recursion formula 

n 

Y £ . x , (f) = + V kf)Y £ . x (f)- (3-7) 

fc=i 

Proof. By the definition (3.4), direct computation leads to 

Ytx,e{f) = / D 1 Y t . x fo,- ,e n {f) + (^l/^-x.fe,- ,o„(f) 

= (V 1 + V 1 f)Y t9Jh ,..,e n (f). 

Similarly, 

Yi-xfiz,- ,e n {f) = {^ > 2 + T^ > 2f)Yt-x,e i ,-,e n {f)- 

Repeating the above arguments then proves the formula (3.7). □ 

Theorem 3. The super Bell polynomials (3.4) possess parity property 

W[(-l) Erj ' +E/lfc We] = (-ir E,) WlW^ (3-8) 
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Proof. From the recursion relation (3.7), it follows that 

n 

W[(-l) Erj+Ew We] = H(-T> k -V k f)Y t . x [{-l)^f r ... 



k=l 



(3.9) 



= (-l) n l[(V k + V k f)Y £ . x [(-l)^f r . x }. 

k=l 

While applying the partitional representation (3.5), we have 

lW(-l) Er 'Yr.*] = {-^YUfr-vi (3.10) 

Hence, combing (3.7), (3.9) and (3.10) proves the formula (3.8). □ 



Theorem 4. The super Bell polynomials (3.4) obey addition property 

1 t\ t m rn 



w/+j/)= E (-ir^--^^ 



(3.11) 



Mi,— ,Mn=0 ri=0 r m =0 i=l 

X ^ / (£-r)-a ; ,(l- M )-0(/)^r. a ;, At .0(5'), 

where n = (1, 2, • • • , n); r[{(l — /x) • n, /x • rc,} \ {0}] denotes the reverse order 
numbers of the n-order permutation {(1 — /x) • n, /x • n}\{0}, which is generated 
from anticommutation of covariant derivatives Vj, j = 1, • • • , n, and obtained 
from a 2n-order permutation {(1 — /Xi)l, • ■ ■ , (1 — /x„)n, /ill, • • • , /x n n} (//j = 
or 1) by taking off all zero terms. The 2n-order permutation is the subscript of 
corresponding covariant derivatives of term y^-rj-aj^i-^-eC/)^-^,^^) kept in 
original order. 

Proof. According to the commutative properties of covariant derivatives (3.1), 
a minus sign in the Leibnitz rule exactly corresponds to an inverse order of the 
n-order permutation {1, 2, • • • , n}. So direct computation shows that 

1 £1 tm rn 

(FG)- i v 1 ...v n di\---d e z(FG)= E (-ir [{(i -^-^ {o}1 E---En( 

M,— ,Mn=0 n=0 r m =0 i=l 

x (FV 1 1 ' f11 ■ ■ ■ V)-^d^- ri ■ ■ ■ d e x 2~ rm F) (GT?f • • • V^d r x \ ■ ■ ■ d r x ™G) , 
which implies (3.11) by replacing F = e* and G = e 9 . □ 
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Let F, G : M™' n — > A be two bosonic super differential functions, then direct 
computation yields 

V 1 V 2 V 3 (FG) = (V 1 V 2 V 3 F)G + {V 2 V Z F)V X G - (V{D 3 F)V 2 G 
+ (V 2 F)(V 1 F)V 2 G + {V{D 2 F)V 3 G - (V 2 F)(ViF)V 3 G 
+ (V 1 F)(V 2 F)V 3 G + F(V 1 V 2 V 3 G), 

in which corresponding six even permutations are {1, 2, 3}, {2, 3, 1}, {3, 1, 2}, {1, 2, 3}, 
{1, 2, 3}, {1, 2, 3}, and two odd permutations are {1, 3, 2}, {2, 1, 3}. 

2.2. Generalized super binary Bell polynomials 

We further define a class of super binary Bell polynomials which play an 
important role in the study of integrability for supersymmetric equations. 

Definition 2. Based on the use of above super Bell polynomials (3.4), the 
super binary Bell polynomials ( ^-polynomials) can be defined as follows 

yi- x ,o(v,w) =Y l . X) e[f r . x ^. e \, (3.12) 

in which we replace the function / and its derivatives by corresponding terms of 
functions w and v respectively, according the following rule 

{Vr-x,re, if EJLi r j + ELi is odd, 
w r . x ^. e , if Ejli r i + ELi is even, 

The super binary Bell polynomials (3.6) is multi-variable polynomials with re- 
spect to various partial derivatives v r . XtfJ ,.g and w r . XtfJ ,.g, rj = 0, • • • ,£j, j = 
0, • • • ,m,fx k = 0, 1, k = 1, • • • ,n. 
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The super binary Bell polynomials also inherit the easily recognizable partial 
structure of the super Bell polynomials. The first few are explicitly calculated as 

y x (v) = v x , y 2x (v, w) = w 2x + vl, y 3x (v, w) = v 3x + Sv x w 2x + v% 
y 01 (v) = T> lV , y 6l e 2 (w, v) = V{D 2 w + (V lV )V 2 v, 
yxfix {v, w) = V x w x + v x V x v, 

(3.13) 

y 2 x,ejy v i w ) = w 2xViv + T>iv 2x + v 2 x T>iv + 2v x V 1 w x , 
y3x,ei(v, w) = v 3x V x v + Zv x w 2x V x v + Zw 2x V x w x + 2>v x V x v 2x 

+ ZvlViWx + V]W 3x . 

We denote the special case of super Bell polynomials by ye. x ,e(v = 0,w) = 
Pi-x,e(w), then it follows from (3.13) that 

P 2x (w) = w 2x , P ix {w) = w Ax + 3wj x , P ei ,e 2 {w) = T>{D 2 w, 

(3.14) 

Pxfii (w) = V x w x , P 3xfil (w) = V x w 3x + ?>w 2x V x w x , 

Theorem 5. The link between super binary Bell polynomials y£. x ,e(v, w) and 
the super Hirota bilinear equation Si • • • S n D x \ • ■ ■ D^F ■ G can be established 
by an identity 

y £ . X!0 (v = InF/G, w = \nFG) = (FG)~ 1 S 1 ■ ■ ■ S n D x \ ■ ■ ■ D^F ■ G. (3.15) 

This formula will be sued to obtain bilinear Backlund transformations of super- 
symmetric equations. 



Proof. Let / = In F, g 
from the definition 2 that 



= In G, then we have v = f — g,w 
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= / + g. it follows 



□ 



y e . X: e(v = lnF/G> = In FG) = Y ( .,.,> U,,po + g s . xM \ 

(3.11) ^_ 1 ^[ { (i_ M ). TliM . n} /0] ^ • • • J] JJ ( * ) 

Hi,— ,Mn=0 n=0 r m =0 «=1 ^ ' 

X F(£-r). a; ,(l- Al ).0[/ S . :C ,A0]^r. 3; , M .0[(-l) ESl+E/tj ^. a; ,A0] 



1 ll lm m / 

(3 = 8) ^ (_l)r[{(l-M).n,/*.n}/0]+E^ir J -+Efc=iM ^ ' ' ' ^ II ( 
/•*!>••• >Mn=0 ri=0 r m =0 t=l 



x F^_ r ). a . j (i_ /1 ). e (/)y T .. a!j/1 . e (^) 

= (FG)- 1 S 1 -~S n L% 1 ..-Djz i F.G. 



For the particular case when F = G, the formula (3.12) reduces to 



G- 2 ^ • • • S n D e x \ ■ ■ ■ DilG ■ G = y t . x , e (0, w = 2\nG) 

0, " + £Jli4 is odd, (3.i 6 ) 



which implies that the P-polynomials can be characterized by an equally rec- 
ognizable even part partitional structure. The formulae (3.15) and (3.16) will 
prove particularly useful in connecting supersymmetric equations with their cor- 
responding super bilinear equations. Once a nonlinear equation is expressible as 
a linear combination of super Bell ^-polynomials or P-polynomials, then it can 
be transformed into a super linear equation. 
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Theorem 6. The super binary Bell polynomials yt. x ,e(v , w) can be separated 
into super P-polynomials and super Bell T-polynomials 

1 h £ m m , x 

y, x Av, W ) = £ ( -i r K(i-M)-n,M-» } /<fl x; • • ■ e n (; ) 

Mi,-,Mn=0 n=0 r m =0 i=l ^ l ' (3.17) 

x P r . x ,n.e(w - v)Y i i_ r) . Xt(1 _ tl) . e (v), 
where only non- vanishing contributions being those for which Y r j + Yl f^k ^ s even 

integer. 

Proof. According Definition 2 of the super Bell polynomials, we have 

y p . x ,u.e(v = 0, w) = 0, as E^' + E Vk is odd ' 

so that by using Theorem 4, 

yt. x ,e( v , w ) = 3W( V > v + q) = Y £ . Xj6 (v + q)\ qp . m , v . g =o 

= E (-ir^-^-^E-Enu) ( 3 - 18 ) 

A*l,---,Mn=0 ri=0 r m =0 i=l ^ ' 

X K-iB,/*.e(9)^-r).as,(l-M)-e( V ) lft,.«, tf .fl=0j 

where q = w — v, the sum + S u k is °dd integer. 
Substituting the relation 

Yr-x,fj,-9\ ( l)\q p . a . :l/ .e=0, Y,Pj+J2 u k is odd = Pr-x,fj,-9\Q) I^Pj+S "k is even ' 

into (3.18) then leads to the formula (3.17). □ 

This theorem implies that the super binary Bell polynomials (3.12) can still 
be "linearized" by means of the Hopf-Cole transformation v — m^, ip — F/G. 
The formulae (3.6) and (3.17) will then provide a way to find the associated Lax 
system of supersymmetric equations. 

Finally, let's through a graph describe general procedure how to use the theory 
of super Bell polynomials that we have developed above. 
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Binlear BT i 



v=\nF/G 



w=lnFG 



SUSY equation : = 

dimensionless field: q 
Bell P(q) system :E(q) = q 21nG > Bilinear form 

Two-fold condition 
E(q) - E(q) = 
constaint 
Binary Bell y(v,w) system 



u=ln ip 



w = hi(r),q), v = h 2 (ri,Q) 
Infinite conservation laws 



Lax pair 



It is clear from this graph to see the close connections among Bell polynomials 
with bilinear equation, bilinear Backlund transformation, Lax pair and conserva- 
tion laws. 



4. The supersymmetric KdV equation 



Consider the supersymmetric KdV equation of Manin-Radul-Mathieu jU E] 

$ t + 3 (®V$) X + = 0, (4.1) 

where $ = t, 9) : -»■ A x is a fermionic super field function with indepen- 
dent variables x, t and Grassmann variable 9. The symbol T> — dg + 9d x denotes 
the super derivative differential operator, which satisfies T> 2 = d x , 9 2 = 0. The 
supersymmetric version of the KdV equation (4.1) describe the time evolution of 
a Grassmann-valued superfield Q(x,t,9) = u(x,t) + 9u(x,t), where u(x,t) is an 
ordinary function and u(x,t) is a Grassmann valued function. The variable x, t 
acquire a Grassmann partner 9, so (x, t, 9) are coordinates in a one dimensional 
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superspace R^' 1 . Since the introduction of the supersymmetric KdV equation 
(1.1) by Manin, Radul and Mathieu [U [6], much attention has been given to its 
mathematical structure and integrable properties. For instances, bi-Hamiltonian 
structure, Painleve property, infinite many symmetries, Darboux transformation, 
Backhand transformation, bilinear form, super soliton solutions and super quasi- 
periodic solutions had been investigated in [TO] [ZE]- Here we see how to apply 
the super polynomials to investigate complete integrability of the supersymmetric 
KdV equation (4.1). 

Theorem 7. Under the transformation $ = 2P(lnG) a ., the supersymmetric 
KdV equation (4.1) can be bilinearized into 

(SD t + SD 3 X )G ■ G = 0. (4.2) 

Proof. The invariance of the equation (4.1) under the scale transformation 

shows that the dimension of the fermionic field $ is —3/2, and it can be related 

2 1 

to a dimensionless bosonic field q : M A ' — > Aq, by setting 

$ = cVq x , (4.3) 

with c G A being free function to be the appropriate choice such that the equation 
(4.1) connects with P-polynomials. Substituting (4.3) into (4.1) and integrating 
with respect to x yields 

E(q) = Vq t + Vq 3x + 3cq 2x V e q x = 0. (4.4) 

Comparing the last two terms of this equation with the formula (3.16) implies 
that we should require c = 1. The equation (4.4) is then cast into a combination 
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form of P-polynomials 

E(q) = P t>e (q)+P 3x>e (q) = 0. (4.5) 

Making a change of dependent variable 

q = 2 In G, $ = 2V(\n G) x 

with G : M 2 ^ 1 — > A , then the property (3.16) shows that the equation (4.5) is 
equivalent to the bilinear equation (4.2). □ 

Starting from the bilinear equation (4.2), it is easy to get super soliton solu- 
tions. For example, the regular one-soliton like solution reads 

$ = £>[ln(l + exp(kx - kH + 0()) x , 

where k G A ,C £ A].. Since solving the equation (4.1) is not our main purpose 
in this paper, the super soliton solutions can be found in details [20J. 

Next, we search for the bilinear Backhand transformation and Lax pair of the 
supersymmetric KdV equation (4.1). 

Theorem 8. Let F be a solution of the equation (4.2), then G satisfying 
(SD X - XS)F • G = 0, 

(4.5) 

(A + L% + 3X 2 D X - 3\D X )F ■ G = 
is another solution of the equation (4.2). This kind of Backlund transformation 
is exactly the same with that given by Liu [TT] 

Proof. Let q = 2 In G, q = 2 In F : IR^' 1 — > A be two different solutions of the 
equation (4.4), respectively, we associate the two-field condition 

E(q) - E(q) = V(q- q) t + V(q - q) 3x + 3q 2x Vq x - 3q 2x Vq x = 0. (4.6) 
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This two-field condition can be regarded as a ansatz for a bilinear Backlund 
transformation and may produce the required transformation under appropriate 
additional constraints. 

To find such constraints, we introduce two new variables 

v — (q — q)/2 — lnF/G, w = (q + q)/2 = InFG, (4.7) 

and rewrite the condition (4.6) into the form 

E{q)-E{q) = 2Vv t + 2Vv 3x + Qv 2x Vw x + 6w 2x Vv x 



= 2V[y t (v) + y 3x {v, w)} + 6R{v, w) = 0, 

with 



(4.8) 



R(v, w) = v 2x Vw x - v x Vw 2x - v 2 x Vv x = Wronskian[^ :Eie (t;, w), y x (v)]. 

In order to decouple the two-field condition (4.8) into a pair of constraints, we 
impose such a constraint which enable us to express R(v, w) as the D-derivative 
of a combination of 3^-polynomials. A possible choice of such constraint may be 

y x ,e(v,w) = \y d (v), (4.9) 

where A G A is an arbitrary parameter. It follows from the identity (4.9) that 

(v x Vv) x = XD e v x - Vw 2x , 

on account which, then R(v,w) can be rewritten in the form 

R(v, w) = X(v x Vv) x - 2Xv x Vv x = X 2 Vv x - XVw 2x - 2Xv x Vv x 

(4.10) 

= v[x 2 y x (v)-xy 2x (v,w)}. 

Then from (4. 7)- (4. 10), we deduce a coupled system of super binary Bell y~ 
polynomials 

y x ,e(v,w)-Xy e (v) = 0, 

(4.11) 

y t (v) + ys x (v, w) + 3X 2 y x (v) - 3xy 2x (v, w ) = o. 



21 

By application of the identity (3.15), under transformation v = In F/G, w = 
In FG, the system (4.11) then leads to the bilinear Backlund transformation (4.5). 
□ 

Theorem 9. The supersymmetric KdV equation (4.1) admits a Lax pair 
{dl + $D - Xd x )p = 0, 

(4.12) 

[Ddt + D8l - 3XVd 2 x + 3(T?$ + X)V + (D<5>)D]p = 0, 
where <p : M 2 ^ 1 — > Ai is a fermionic eigenfunction. 

Proof. By transformation v = lnip, it follows from the formulae (3.6) and 

(3.17) that 

y e (v) = Vrj}/^, y x , e (v, w) = Vq x + 

ytiv) = ^tM y2x(v,w) = q 2x + ^2x/^: ^3x(^, w) = ^2x^x1^ + ^x1^, 

on account of which, the system (4.12) is then linearized into a Lax pair with a 
parameter A 

L\ip = (Dd x -XV + Vq x )tP = 0, 

L 2 iP = (d t + dl + 3q 2x d x + 3X 2 d x - 3Xd 2 + q 2x )^ = 0, 
which is equivalent to the formula (4.12) b by replacing Vq x with $, and tp with 

Dip. It is easy to check that the integrability condition of the Lax pair 

[L u L 2 }^ = 

is satisfied if $ is a solution of the supersymmetric KdV equation (4.1). □ 

Finally, we show how to derive the infinite conservation laws for super KdV 
equation (4.1) based on the use of the binary Bell polynomials. 

Theorem 10. The supersymmetric KdV equation (4.1) possesses the follow- 
ing infinite conservation laws 

I n , t + DF n = 0, n = l,2,-. •. (4.13) 
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where the fermionic conserved densities I' n s are explicitly given by recursion re- 
lations 



h = Vq x = $, I 2 = I l)X = $ x , 

n 

I n +i = I n ,x + ^ hT^In-k, n = 2,3, 



(4.14) 



fc=i 



and the bosonic fluxes F^s are given by recursion formulas 
F 1 = V<$> 2x + 3$$, + 3(P$) 2 , 

F 2 = X?$ 3;E + 3($$ 2a; + + 6V$V<f> x: 

n 

F n = VI n , 2x + 3 Y^ihln+i-k + VI k VI n+1 _ Kx ) + 3P$D/„ 



(4.15) 



fc=i 

+ VhViph, n = 3,4,-... 

i+jr'+fc=n 

Proof. The conservation laws actually have been hinted in the two-filed 
constraint system (4.9)-(4.11), which can be rewritten in the conserved form 

y xfi (v,w)->y e {v) = Q, 

(4.16) 

d t y e {v) + V[y 3x (v, w) + 3A 2 34(t;) - sxy 2x (v, w)] = 0. 
by applying the relation T>y t {y) = d t yeiy) = Vv t . 

By introducing a new fermionic potential function 
V =(Vq-Vq)/2: R 2 / -+ A u 
it follows from the relation (4.8) that 

Vv = r], Vw = r] + Vq. (4.17) 
Substituting (4.17) into (4.16), we get a super Riccati-type equation 

rj x + rjVrj + Vq x -\r] = 0, (4. 18) 

and a divergence-type equation 

T] t + V[Vr) 2x + 3\r]r]x + Sq 2x Vr] + 3Vr]Vr] x + (Vr]) 3 } = 0, (4.19) 
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where we have used the equation (4.18) to get the equation (4.19). 
To proceed, inserting the expansion 

oo 

r / = ^/ n (Pg,g x ,---)A- ra , (4.20) 

71=1 

into the equation (4.18) and equating the coefficients for power of A, we then 
obtain the formulas (4.13). 

Finally, substituting (4.20) into (4.19) yields 



n=l 



oo 



n=l n=l n=l n=l 



0. 



oo oo oo 

\3 



+3 J2 VI n^ E VI n,^ H + E VI n^f 
n=l n=l n=l 

which leads to infinite consequence of conservation law equation (4.13) by equat- 
ing the coefficients for power of A. □ 

It follows from the conservation equation (4.13) by using (2.2) that 



J j I n dxd9^j =~JJ (DF n )dxd0 = 0, 



which implies that I' n s are fermionic conserved densities. We present recursion 
formulas for generating an infinite sequence of conservation laws for each equa- 
tion, the first few conserved density and associated flux are explicit. The first 
equation of conservation law equation (4.13) is exactly the supersymmetric KdV 
equation (4.1). In conclusion, the supersymmetric KdV (4.1) is complete inte- 
grable in the sense that it admits bilinear Backlund transformation, Lax pair and 
infinite conservation laws. 

5. The supersymmetric sine-Gordon equation 
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The classical sine-Gordon equation 

(p xt = sin(j) (5.1) 

has applications in various areas of physics including nonlinear field theory, solid- 
state physics, nonlinear optics, elementary particle theory and fluid dynamics, see 
[38]- [12] and references therein. The supersymmetric extension of the equation 
(5.1), i.e. the supersymmetric sine-Gordon equation [43j-[50j 

V 1 V 2 $ = sin$ (5.2) 

is constructed on the four dimensional superspace (x, t,e x ,e 2 ) e RZ . Here, 
$ = $(x, t, 6\, 62) '■ M^' 2 — > Ao is a scalar bosonic superfield; The variables x and 
t represent the even coordinates on the two-dimensional super- Minkowski space, 
while the quantities 9\ and 82 are anticommuting odd coordinates which satisfy 
the anticommutation relations 

91 = 91 = 0, {9 1 ,9 2 }=0. 

The T>\ = de 1 + 9\d x and T>2 = de 2 + #2<9t are two covariant derivatives which 
satisfy the anticommutation relations 

v\ = d x , v 2 2 = d u [v u v 2 ] = o. 

The supersymmetric version of the sine-Gordon equation was introduced from 
purely physical motivations [13] . It is becoming increasingly interesting to inves- 
tigate the supersymmetric sine-Gordon equation because of its close relation to 
string theories and statistical physics [44]- [16]. In recent publications, a super- 
space extension of the Lagrangian formulation has been established for the super- 
symmetric sine-Gordon equation [18]. The bilinear method is used to construct 
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multi-super soliton solutions jl?]. The supersymmetric sine-Gordon equation ad- 
mits a Lax pair, and a connection was established between its super-Backlund 
and super-Darboux transformations [HJH9]. The method of symmetry reduction 
is systematically applied in order to derive invariant solutions of the supersym- 
metric sine-Gordon equation [50]. The prolongation method of Wahlquist and 
Estabrook was used to find an infinite-dimensional superalgebra and the associ- 
ated super Lax pairs [5T] . 

Here we study the integrable properties of the supersymmetric sine-Gordon 
based on the use of generalized super Bell polynomials. The bilinear form, bilinear 
Backlund transformation, Lax pair and infinite conservation laws systematically 
are obtained with our method. 

Theorem 11. Under the transformation 

$ = 2i\n(F/G), 

the supersymmetric sine-Gordon equation (5.2) admits the bilinear form 

2S 1 S 2 F ■ F + G 2 = 0, 2S 1 S 2 G ■ G + F 2 = 0, (5.3) 

where F, G : R][ A are two bosonic functions. 

Proof. As before, the invariance of the supersymmetric sine-Gordon equation 
(5.2) under the scale transformation 

x ^ Ax, t A _1 t, 0! -> A 1/2 #i, 6 2 -> \~ 1/2 6 2 , $ ->■ $ 

shows that the dimension of the bosonic superfield $ is zero, and so we may 
introduce a dimensionless bosonic field q by setting 

$ = cq, (5.4) 
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in which c G A is free constant to be determined. Substituting (5.4) into (5.2) 
yields 

2V l V 2 q = P ei02 (p + q) - P ei e 2 (p - q) = i{e~^ - e ic <*)/c (5.5) 

where p : M^' 2 — > A is an auxiliary function. If one chooses the constant c — 2i, 
the equation (5.5) is then cast into a linear combination form of P-polynomials 

2iWp + q)~ 2P 9l e 2 (P ~ q) + exp(-2g) - exp(2g) = 0, 

which can be decoupled into a system 

E^p, q) = 2P eid2 (p + q)+ exp(-2g) = 0, 

(5.6) 

E 2 (p, q) = 2P dl e 2 (p - q) + exp(2g) = 0. 
Multiplying the first equation by exp(p + q), the second equation by exp(p — q) 

in the equation (5.6) yields 

2 exp(p + q)Pe 1 9 2 (p + q) + exp(p - q) = 0, 

(5.7) 

2 exp(p - q)P eid2 (p - q) + exp(p + q) = 0. 
By transformation 

q = \n(F/G), p = ln(FG) <^ $ = 2iq = 2i \n(F/G), p = \n(FG) 

and using the property (3.16), then the equation (5.7) gives the bilinear form 
(5.3) for the supersymmetric sine-Gordon equation (5.2). □ 

Theorem 12. Let (F,G) be a solution of the equation (5.3), then (F,G) 
satisfying 

S 1 G-G = XgFF, S ± F ■ F = -XgGG, 

S 2 F.G = -L 9 GF, S 2 G.F = -±. 9 FG, (5g) 
„ / FF GG\ „ 1 (GF FG 
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is another solution of the equation (5.3), where g : M^' 2 — > Ai fermionic auxiliary 

superfield and A G A is even parameter. 

Proof. In order to obtain the bilinear Backlund transformation and Lax pairs 

of the equation (5.2), let p,q and p, q be two solutions of the equation (5.6) and 

consider the associated two-field condition 

E ± (p, q) - E 1 {p, q) = 2V 1 V 2 (p - p) - 2V{D 2 {q- q) 



+ e q+q {e q - q - e q - q ) = 0, 
E 2 (p, q) - E 2 (p, q) = 2V 1 V 2 (p - p) + 2V 1 V 2 (q - q) 

+ e ~ {q+q \e q - q -e q - q ) = 0, 

where 

p = \n(FG), q = \n(F/G), 

We introduce variables 

Vl = HG/G), v 2 = \n(F/F), v 3 = \n(F/G), v A = \n(G/F), 

Wl = \n(GG), w 2 = ln(FF), w 3 = ln(FG), w 4 = ln(GF), 
from which, we have relations 

q - q = v 2 - v 1 = w 3 - u> 4 , q + q = v 3 - v 4 = w 2 - Wi, 

p - p = Vi + v 2 = v 3 + v 4) p + p = Wi+w 2 =w 3 + w 4 

and 



(5.9) 



(5.10) 



v 1 =v 4 + q, v 2 = v 3 -q, w 1 = w 4 -q, w 2 =w 3 + q. (5.11) 

By using the mixed variables (5.10), it follows that from (5.9) 
AV{D 2 v x + e V3 - V4 (e V2 ~ vl - e Vl ~ V2 ) = 0, 

(5.12) 

AV{D 2 v 2 + e VA ~ Vz (e Vl - V2 - e V2 - Vl ) = 0, 
which may produce the required bilinear Backlund transformation under an ap- 
propriate additional constraint. We choose a constraint 

V 1 v 1 =y ei (v 1 ) = Xge V3 ~ v \ (5.13) 
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where g : M 2 ^ 2 — > Ai fermionic auxiliary superfield and A G A is even parameter. 
The fermionic function g is introduced because of supersymmetry and the odd- 
ness of the superspace derivatives T> 1: T> 2 . The constraint (5.13) reduces the first 
equation in (5.12) into 

-AXV 2 g - 4XgV 2 (v 3 - v 4 ) + e V2 ~ Vl - e Vl ~ V2 = 0. (5.14) 

Since the term V 2 (v 3 — t> 4 ) should be fermionic function, we make a constraint 

V 2 (v 3 -v 4 ) = ^gh, (5.15) 

where h is a bosonic function to be determined. On account of this constraint, it 
follows from (5.14) that 

Vi9 = ^(^-e^ 2 ), (5.16) 

which holds because g 2 = 0, g being fermionic. 

By means of the system (5.15) and (5.16), the second equation in (5.12) reads 

V 2 (V lV2 + \ge v *- vs )=0, 

which is satisfied if we choose 

y 01 (v 2 )=V 1 v 2 = -\ge v *- V3 . (5.17) 

On the one hand, using the relation (5.11), we have 
AViV^Vs - Vi) + 4£> 2 £> 1 (t; 1 - v 2 ) = 8V ± V 2 q 

(5.18) 

_ _ gUl— ^^gUS— U4 _|_ e V4-V3^ _|_ _|_ gUl —V2 ^gl^-lM _ gl>4-U3^ 

On the other hand, it follows from (5.13), (5.16) and (5.17) that 
AV 2 V,( Vl - v 2 ) = 4\(V 2 g)(e v *- v * + e") 

(5.19) 



Combining (5.15), (5.18) and (5.19) yields 



AV 2 V 1 (y 3 - v 4 ) = \{V ig )h = (e v *- Vi - e v *- V3 )(e V2 - Vl + e^ 2 ), 
A 



which implies that we may choose 
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V l9 = \(e V3 - V4 - e V4 ~ V3 ) 



(5.20) 



and 



Thus, we have 



h = (e V2 ~ Vl +e Vl ~ V2 ). 



T> 2 (v 3 -v 4 ) = —g{< 



V2—V1 I „V1— V2 



which can be written as a pair of ^-polynomials 



^2(^3) = V 2 v 3 = -^ge 



(5.21) 



Combining (5.13), (5.16), (5.17), (5.20) and (5.21) gives bilinear Backlund 
transformation (5.8) of the supersymmetric sine-Gordon equation. □ 

Finally we derive Lax pair of the supersymmetric sine-Gordon equation. 
Theorem 13. The supersymmetric sine-Gordon equation (5.2) admits a Lax 

pair 



= Mm 



V 2 ^> = 



1 <i> \g X 

1 1 

y 2 / 



-Ix 9e 



7 



together with 



Vi9 = *[-, r ' = T e — - e — 

>4 ^3/ A V ^4 ^3 



(5.22) 



where * = (^3, ^4) 1 
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Making use of the Hopf-Cole transformation 

v 3 = \nip 3 , w 4 = ln^ 4 , 

then the system (5.13), (5.16), (5.17), (5.20) and (5.21) can be linearized into a 
Lax pair (5.22). It is easy to check that the integrability condition 

V 2 M + VxN - [M, N}=0 

is satisfied if $ is a solution of the sine-Gordon equation (5.2). 
If we choose a transformation 

0i = fa =ipl, 9= .f >3 ; 

2ifaip 4 

then the Lax pair (5.22) is also equivalent to a linear system in 3 x 3 matrix form 



1 



'4Z>i$ A 



V<$1 = - -4X>!$ -A | n, 
1 1 -4A 4A 




(5.23) 



where f2 = (fa, 2 , 03) T , (fix, fa '■ K/{ 2 — > A are bosonic functions and fa : R^' 2 — > 
Ax is a fermionic function. The system (5.23) also can be obtained from (5.13), 
(5.16), (5.17), (5.20) and (5.21) by setting 

01 n/ \ r>-x. , i 01 03 



2(v 4 - U3 ) = ln^i, 2(u 2 -« 1 ) = -2z$ + ln-ji, 3 



02 ' " 02 2ly/fafa 

The compatibility of the linear system (5.23) in superspace is equivalent to the 
equation (5.2). The system (5.23) is the same as obtained in [UJ, but here it is 
derived systematically from the super Bell polynomials and Lax pairs. □ 
Noting the transformation relation 

Vl -v 2 = i($ - $)/2, v 3 -v i = + $)/2, 
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then it follows from equations (5.13), (5.17), (5.18), (5.21) and (5.22) that 



£>i($ - <3>) = Agcos 



$ + $ 
~2~ 



P 2 ($ + $) = i^cos^-^), (5.24) 

, . . $ + $\ ^ 1 . 
Pi5 = Asm | — ^— I , £> 2 = ^ sin I — 



which is the Backlund transformation of the supersymmetric sine-Gordon equa- 
tion. The compatibility of the Backlund transformation (5.24) is the supersym- 
metric sine- Gordon equation for both $ and <E> separately. The super Backlund 
transformation (5.24) reduces to the classical Backlund transformation of the 
purely bosonic sine-Gordon equation when fermions are equal to zero. 



6. Concluding Remarks 



In this paper, we have introduced a class of super Bell polynomials which play 
an important role in the characterization of bilinear Backlund transformation, 
Lax pairs and infinite conservation laws of supersymmetric equations. To the 
knowledge of the authors, this is the first work on the super Bell polynomials and 
their applications to super integrable systems. We believe that there are still many 
interesting deep relations between generalized Bell polynomials and integrable 
structures, which remain open and worth to be considered. For instance, (i) How 
to explore the relations between the super Bell polynomials with symmetries, 
Hamiltonian functions, etc. (ii) How to define a class of discrete Bell polynomials 
and apply them in discrete equations. We have some ideas on these questions 
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and will intend to return to them in some future publications. 
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